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I. INTRODUCTION 

In a recent paper Ref. [lj, hereafter referred to as I, we 
considered how the new plateaus in the Hall conductivity 
of graphene discovered by Y. Zhang et al£ develop, as- 
suming the weak-coupling regime of the magnetic cataly- 
sis scenario. We have become aware of the experiments 3 
showing that the v — plateau in the Hall conductivity 
exists even at rather high temperatures. This fact in- 
dicates the relevance of a strong-coupling regime. It is 
the purpose of the present work to complete the theoret- 
ical analysis of I by considering also the strong-coupling 
regime. 

A graphite monolayer, or graphene, has become a new 
exciting topic in physics of two dimensional electronic 
systems.— A qualitatively new feature of graphene is 
that the eigenfunctions and the eigenvalues of the low 
energy quasiparticle excitations are described by the rel- 
ativistic 2+1 dimensional Dirac theory. The spinor struc- 
ture of the wave functions is a general consequence of the 
honeycomb lattice structure of graphene with two carbon 
atoms per unit cel\2& When a magnetic field is applied, 
noninteracting Dirac quasiparticles occupy Landau levels 
(LLs) with the energies 

E n = sgn(n)^2\n\hv 2 F \eB\/c 
w 424 sgn(n) a/R ^B [T] K, n = 0, ±1, ±2, . . . 

Here to estimate the energies E n , the value of the Fermi 
velocity vf — 10 6 m/s was used, and the magnetic field 
B orthogonal to the graphene's plane is given in Tesla. 
Several anomalous properties of graphene are attributed 



to the presence in the spectrum (flj of the n = field 
independent lowest Landau level (LLL). In the diagonal 
conductivity, the anomaly manifests itself as the phase 
shift of 7r of the quantum magnetic oscillations expected 
theoretically both using the semiclassical quantization 
condition for the quasiparticles with a linear dispersion 9 
and from a microscopical calculation valid both for the 
massless and massive Dirac fermions i 10 : 11 Accordingly, in 
the Hall conductivity the anomaly results in the anoma- 
lous integer quantum Hall (QH) effect with the plateaus 
at the filling factors v = ±4(|n| + 1/2)^^^^ These 
theoretical conceptions allowed to identify unambigu- 
ously the Dirac quasiparticles in the two independent 
experiments! 16 ' 17 At present a lot of studies on graphene 
are concentrated on phenomena whose understanding de- 
mands going beyond unconventional but yet rather sim- 
ple physics of noninteracting Dirac quasiparticles. 

On the experimental side, for a magnetic field B > 
20 T, the appearance of additional QH plateaus with 
the filling factors v — 0, ±1,±4 was reported in Ref. 0. 
The theoretical studies of these additional plateaus in 
graphene can be divided into four classes, 
(i) Fractional QH effect. Although there is no experi- 
mental evidence for such an effect in graphene so far, 
there are a few theoretical papers, where this possibility 
is discussed^ii&i 9 -^ 

The remaining three scenario consider various pos- 
sibilities of breaking the U(4) symmetry of the non- 
interacting Hamiltonian of graphene. While this sym- 
metry remains intact in the presence of a long-range 
Coulomb interaction and at a nonzero chemical potential 
fi, it is explicitly broken in the presence of the Zeeman 
and some types of short-range interactions. 
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(ii) A model with local (on-site) interactions which ex- 
plicitly break the U(A) symmetry was considered in 
Ref. E|. 

(iii) An analogy between the four-fold degeneracy of LLs 
in graphene associated with the U(4) symmetry and the 
SU(4) ferromagnetism studied previously in the bilayer 
quantum Hall systems^ is exploited in Refs. I23II24I In 
this QH ferromagnetism scenario the QH plateaus with 
all integer values of the filling factor v occur. The current 
experimental data, 2 however, do not seem to support the 
existence of the plateaus with v = ±3, ±5, . . .. 

(iv) The magnetic catalysis scenario was considered in 
I and in Refs. I25II26I . It is based on the phenomenon 
of the electron-hole (fermion-antifermion) pairing in a 
magnetic field revealed in field theory 2 ^ and the anal- 
ysis in Refs. [28. 29 30 3l|, where this phenomenon was 
considered in graphene. (This analysis was originally in- 
spired by the early experiments in hig hly oriented py- 
rolytic graphite. 32 ) In Refs. 28. 29. 30. 3_l|, the spontaneous 
breakdown of the U(4) symmetry by generating a dynam- 
ical excitonic gap A was considered. A new development 
for this scenario suggested in I was to fit the data of 
Ref. by including both the Zeeman term and the exci- 
tonic gap. The central feature of this scenario is that the 
only plateaus in the Hall conductivity a xy are those with 
v = 0, ±1 and v — ±2k (k = 1,2,.. .), i.e., the plateaus 
observed in experiment. 2 . While in I the excitonic gap 
is produced by the Coulomb interactions, local (on-site) 
interactions are used for this purpose in Ref. [25J. A dy- 
namics relating to the magnetic catalysis scenario was 
considered in Ref. [HI, in which an excitoniclike gap is 
produced by electron-phonon interactions. 

As was emphasized in I, while the plateau v — could 
even appear either due to the spin splitting or the exci- 
tonic gap A alone, the plateaus v = ±1 arise only if both 
the spin splitting and the gap A are non-vanishing. In 
other words, the plateaus v — ±1 are generated by the 
dynamics which completely removes the U(4) degeneracy 
of the LLL. This can for example be seen explicitly from 
the expression for the Hall conductivity due to the n = 
Landau level in the clean limit (i.e., in the limit of the 
vanishing scattering rate of quasiparticles): 

e 2 

v X y = — ^-sgn(eB) ^ 
x [sgn( M+ )0(|/i + | - A) + S gn(/O0(| M _| - A)] , 

where (i± = fi ± Ez with Ez being the Zeeman energy. 
The fitting procedure in I is heavily based on the assump- 
tion that the Zeeman energy is 

E z = ^» B B~0.67B[T}K, (3) 

where [Ib = eh/(2mc) is the Bohr magneton and the 
Lande factor in graphene is ql — 2. For typical strengths 
of the magnetic field used in the experiment, B < 45 T, 
the Zeeman energy © is Ez < 30 K. By combining this 
observation with the fact that the v = and v = ±1 



plateaus have comparable widths^ it was concluded that 
the excitonic gap A is of the same order as Ez- Indeed, 
as seen from Eq. ((2|), it is the interplay between these 
two energy scales, A and Ez, that determines the size of 
the lowest plateaus. The required magnitude (of several 
dozens of Kelvin) for the excitonic gap corresponds to 
the weak-coupling regime. In J, the best fit was found at 
g = e 2 /(ehvp) ~ 0.07, where e is the dielectric constant of 
the medium. An immediate implication of such a weakly 
coupled dynamics is that the typical magnitude of the 
critical temperature, at which the gap disappears, is of 
the same order of a dozen of Kelvin. 

In part, the motivation for the present work is the new 
experimental data 3 - suggesting that the v = plateau, 
observed for strong magnetic fields, may persist even 
at rather high temperatures. This alone would suggest 
the relevance of the strong-coupling regime with 5 > 1. 
Moreover, as will be shown in Sec.|TT]in detail, this regime 
is also consistent with the structure of the higher plateaus 
in the Hall conductivity reported in Ref. 0. 

In contrast to the weak-coupling regime, the dynam- 
ics at strong coupling is non-perturbative and, therefore, 
permits no rigorous quantitative analysis. In addition, 
the simple type of the dynamics discussed in I is not 
sufficient to fit the experimental data in Ref. d. This 
is because of a very large hierarchy between the energy 
scales set by the Zeeman splitting © and the excitonic 
gap at strong coupling. In a consistent approach, as will 
be discussed below, an anomalous enhancement of the 
Zeeman-like splitting will be required. 

Because of the difficulties due to the non-perturbative 
dynamics at strong coupling, we use a phenomenological 
approach in this paper. Therefore, most conclusions of 
the analysis below should be viewed only as qualitative. 
As we shall see, however, many of them are likely to be 
very robust. For example, the critical temperature for the 
excitonic condensate as well as for the related appearance 
of the plateaus v = 0, ±1 should be of the same order as 
the excitonic gap, i.e., a few hundred Kelvin. 

The paper is organized as follows. In Sec. [Til we ana- 
lyze the experimental data of Ref. |2| in detail, and extract 
the constraints on the dynamics to be implemented in a 
phenomenological model. The model itself is introduced 
in Sec. IIII1 In Sec. IIV1 we present a fit of the experimen- 
tal data for the v = 0, ±1 as well as the higher plateaus, 
and underline the specific features of the strong-coupling 
regime. In Sec. [Vj for completeness and better under- 
standing the role of the Zeeman splitting, we consider the 
dynamics in the strong-coupling regime with no enhanced 
Zeeman splitting. In this case, besides the standard se- 
quence v = {An + 2), the only additional plateau is that 
with v = 0. In Sec. IVII the main results of the paper are 
summarized. 
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II. ANALYSIS OF THE EXPERIMENTAL DATA 

The low-energy quasiparticlcs excitations in graphene 
are described in terms of a four-component Dirac spinor 
= {^KA<7,ipKB^,^K'Ba,ipK'Aa)- This spinor com- 
bines the Bloch states with spin a = ±1 on the two dif- 
ferent sublattices (A, B) of the hexagonal graphene lat- 
tice and with momenta near the two inequivalent points 
(K, K') at the opposite corners of the two-dimensional 
Brillouin zone. 

The free, low-energy quasiparticlc Hamiltonian for zero 
carrier density (or \x = 0) and in absence of the Zeeman 
splitting can be recast in a relativistic form, 

Hq = -iv F J d 2 r^ a (7 1 fiV x + l 2 hW y ) (4) 

where ^ a — ^^j is the Dirac conjugated spinor and 
summation over spin a is understood. Notice that the 
Fermi velocity vf ~ 10 6 m/s plays the role of the speed 
of light. In Eq. ([4]), 7", v = 0, 1,2, are 4x4 gamma 
matrices belonging to a reducible representation of the 
Dirac algebra: 7^ = T3 ® (T3, 1T2, —iT%), where the Pauli 
matrices f, r act in the subspaces of the valley (K,K') 
and sublattices (A, B) indices, respectively. The matrices 
satisfy the usual anticommutation relations {7^,7"} = 
25^, = = 0,1,2. The orbital ef- 

fect of a magnetic field B applied perpendicular to the 
graphene plane is included via the covariant derivative 
V = d+ (ie/He)A. 

The explicit form of the interaction Hamiltonian and 
the effective action are derived and discussed in I. We de- 
scribe their features in the strong-coupling regime below. 

A. Relationship between the gate voltage and p 

In experiment s] 2 ' 4 ' 16 ! 17 the density of carriers (or carrier 
imbalance p which is the difference between the densities 
of electron and holes) is tunable by the gate voltage V g 
applied to the Si substrate of a graphene device. The 
measurements of the Hall coefficient and a (consistent 
with them) theoretical estimate give the relationship 16 

p = a{V g -V ), a « 7.3 x 10 10 citT 2 V~\ (5) 

where the shift Vq is attributed to the shift of the Dirac 
point due to an unintentional doping* 4 - In our analysis of 
the experimental data^ the value of the shift is in the 
range from 0.8 V to 5.8 V. 

On the theory side, the variation of the density of car- 
riers in the system is modeled through adding to the 
Hamiltonian (g]) the term -^^7°*^ = -p$>\ with a 
tunable chemical potential p. In order to make a connec- 
tion between the experiment and the theory, one needs 
to know the relationship between the chemical potential 
\x and the carrier imbalance p. 

We have examined the data of Ref. d relying on the 
following two natural assumptions: 



(i) the middle point of a step between two neighboring 
plateaus is associated with the value of the chemical 
potential in the middle of a smeared Landau level, 
and 

(ii) the middle of a plateau is associated with the value 
of the chemical potential in the middle of the energy 
gap between two Landau (sub)-levels. 

The meaning of these two assumptions might be easier 
to grasp by referring to the correspondence between the 
quasiparticle energy spectrum and the structure of the 
Hall conductivity plateaus as shown, e.g., in Fig.[5]below. 

We note that the structure of the higher plateaus 
can be used to extract the dependence p.(V g ). In the 
corresponding analysis it is useful to assume that only 
the plateaus v — 0,±1,±2 (connected with the LLL) 
are affected by the excitonic gap. This qualitative fea- 
ture is a rigorous outcome in the weak-coupling regime 
of I. If it holds also at strong coupling, which is at 
least plausible^— the dynamics of the exciton conden- 
sation would not contaminate the structure of the higher 
plateaus in the Hall conductivity. Then, for the filling 
factors v = ±3, ±4, ±5, one derives the following rela- 
tions: 

V g (v = ±A) = V g (fi = E ±1 ), (6a) 
V g (y = ±3) = V g (p = E ±1 t E z ), (6b) 
V g (v = ±5) = V g (p = E ±1 ± E z ), (6c) 

where Eqs. (|6b|) , (f6"c|) , and Eq. (|6"aj) are based on the as- 
sumptions (i) and (ii), respectively [see also Fig.[5]below]. 
In the analysis it is reasonable to accept that the Zeeman- 
like energy Ez, characterizing the Landau level splitting, 
is considerably smaller than E\ (i.e., Ez even if 

Ez is substantially enhanced compared to the ordinary 
Zeeman energy in Eq. ([3]). As we will see in Subsec. IIB 
and IIC below, this hierarchy of scales is indeed consis- 
tent with the strong-coupling regime: while Ez — 100 K, 
the energy E\ ~ 1000 K. The results of the analysis are 
summarized in Figs. [T] and [2] 

In Fig. [TJ the experimental data for three specific com- 
binations of the gate voltage differences are plotted as 
functions of the external magnetic field. After taking 
into account the correspondence between the filling fac- 
tors and the values of the chemical potential in Eq. ©, 
we see that all three gate voltage differences, 
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determine the average value of the voltage (measured 
from Vq) that, up to higher order corrections, corresponds 
to p, = E x , i.e., fj. = 424^5 [T] K, see Eq. (TI|). (Note 
that to linear order the corrections due to Zeeman split- 
ting cancel in all three combinations.) From Fig. [TJ we 
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FIG. 1: (Color online) The compilation of the experimen- 
tal data from Ref. for the gate voltage V g {v) differences, 
l/2{V g {u = 4) - V g {y = -4)], l/2[V B {u = 5) - V B (y = -3)] 
and l/2[V g (u = 3) — V g {v = —5)] as a function of the magnetic 
field B. 



see that the dependence of the voltage on the magnetic 
field is described almost perfectly by the following lin- 
ear fit: V g (B) « 1.38 B[T] V where V g = \V g - V \ and 
Vq describes the shift of the Dirac point, see Eq. ([5]). 
Therefore, by trading the magnetic field B for the corre- 
sponding chemical potential, i.e., u = 424-^/ B [T] K, we 
get 

By inverting this relation, we arrive at the following one: 
fi = KSgn(V g -Vo)\/\V g -V \, K w361KV- 1/2 . (9) 

The relations (8J and ([9]) will be used below to make a 
contact between the theory and the experiment. [The 
relation @ is shown in Fig. [5] below by a dotted line.] 

Another independent derivation of the relationship be- 
tween V g and [i can be obtained by fitting a sequence of 
several Landau levels E n at a fixed value of the mag- 
netic field. In Ref. d, several lowest plateaus that cor- 
respond to the filling factors v = An (with an inte- 
ger n) were observed for the four lowest values of the 
field. In accordance with our assumption (ii), the cor- 
responding values of the voltage are associated with the 
chemical potentials in the middle of the energy gap be- 
tween the two sublevels of the Zeeman-split nth Landau 
level, i.e., /i = E n to linear order. By making use of 
Eq. ([T]), the corresponding chemical potentials are given 
by a = 212sgn(i/) v /|j/B| [T] K with v = An. 

The experimental data for the voltage differences that 
correspond to the first few observed sets of the filling 
factors are plotted versus vB in Fig. [2j This depen- 
dence is well approximated by a linear function, V g ~ 
(1.05 + 0.33 vB [T]) V. Note that because of limited data 
for low B, the value 1.05 V of a small intercept term 
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FIG. 2: (Color online) The compilation of the exper- 
imental data from Ref. @ for the voltage differences 
l/2[V g (^/\v/2\hv F \eB\/c) - Vg(- y / \v/2\liv F \eB\/c)] that cor- 
respond to the first few observed sets of the filling fac- 
tors as a function of vB for four values of the field B — 
9,11.5, 17.5, 25T. 



< Vo is not reliable and will be omitted. Then, by ex- 
pressing vB in terms of the chemical potential (as dis- 
cussed in the previous paragraph), we arrive at the re- 
lation V g ~ (/i[K]/369) 2 V, which agrees reasonably well 
with Eq. ([8]). 

Interestingly, the relation §§§ with practically the same 
coefficient k also follows from the equation a 2 = nh 2 v F \p\ 
for the ideal gas of Dirac quasiparticles at T = B = 
used in Refs. Il4ll29l . Note however that the simultane- 
ous use of Eq. ([9]) and the usual (i.e., without enhance- 
ment) Zeeman splitting ([3]) would lead to inconsistencies 
in fitting the experimental data. This was the reason for 
considering another relation between fi and V g at weak- 
coupling in I. 

B. Zeeman term 

Now we are in a position to discuss the second, Zee- 
man term hbB^j ^ 1 ^ = ubB^ which has to be 
added to the free Hamiltonian (|4]) (here the 03 matrix 
acts on spin indices). Considering the weak-coupling 
regime in I, we assumed that the Zeeman splitting is 
usual as given by Eq. ([3]). There are, however, theoretical 
argument a 23 ! 24 ' 34 that the Coulomb interaction in the ex- 
change channel may strongly enhance Zeeman splitting 
and lead to a spin gap Az (expressed through the con- 
densate (^7°o"3^ r )) as large as a few hundred Kelvin. 

Using the experimental data of Ref.0, the energy of the 
Zeeman splitting can be found from the size of v = ±4 
plateau whose appearance is attributed to the lifting of 
the spin degeneracy of the n = ±1 Landau level. We es- 
timate the size of this plateau by extracting the following 
gate voltage differences: SV g = l/2[V g (p = E\ + Ez) — 
Vg(ji = E 1 - Ez)] w l/2[V g {a = £_j + E z ) - V g (p = 
E-i — Ez)]- The results are presented in Fig.PJJ We find 
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FIG. 3: (Color online) The compilation of the experimental 
data from Ref. @ for the size of the v = 4 plateau that deter- 
mines the Zeeman splitting of the first Landau level. 



that the size of a xy — ±4e 2 /h plateau is reasonably well 
fitted by the dependence 

6V g (B) m (-1.1V-B[T] +0.41B[T])V (10) 

for 9T < B < 45 T. On the other hand, if dpi, cor- 
responding to SV g , is reasonably small, one can dif- 
ferentiate expression ([8]) to obtain the relation SV g « 
2^[K]<5/i[K]/(361) 2 V (here the differentials dV g and dpi 
were replaced by 5V g and 8 pi, respectively). Now, sub- 
stituting the chemical potential pi = E\ = A2Ay/B[T] K 
and the splitting dpi — Ez in this expression for SV g and 
comparing it with Eq. (|10p . we arrive at the relation 

E z = (63-/B[T] - 169) K. (11) 

This is 3 to 8 times larger than the usual Zeeman energy 
([3]) for 9 T < B < 45 T, providing another argument in 
favor of strong coupling. The large difference between the 
enhanced Zeeman energy (|11[) and conventional Zeeman 
energy ([3]) is likely caused by dynamical effectsi 23 ' 24 ' 34 In 
general, these effects should be sensitive to the tempera- 
ture and chemical potential, so that Ez is also a function 
of T and pi. Here, however, we restrict ourselves by the 
dependence Ez on B. 

C. Gap terms and the v = ±1 plateaus 

As is clear from Eq. @, the v = ±1 plateaus of the 
Hall conductivity observed in Ref. can exist only if there 
is Zeeman splitting and the spectrum ([!} is gapped, viz. 

E n = sgn(n)^2\n\hv 2 F \eB\/c+A 2 . (12) 

In this case the four-fold degeneracy of the LLL (n = 0) 
is completely lifted. 

Considering the dependence of a xy on pi in Eq. one 
can see that while the size of the v = ±1 plateaus is 5pi\ ~ 



2Ez, the size of the 0th plateau is Spio ~ 2(A — Ez). The 
experiment in Ref. [U measures the dependence a xy (V g ) 
and indicates that the sizes of the 0th and v = ±1 
plateaus (as a function of V g ) are comparable. Using 
Eq. ([9]), we estimate the V^-size of the v = ±1 plateaus 
as AAEz I k 2 (in Volts), while the size of the 0th plateau is 
2(A — Ez) 2 /k 2 (in Volts). By simply requiring these two 
quantities to be equal, we derive A ~ (2 + y/3)Ez- By 
making use of the Zeeman energy (fTTj) extracted from ex- 
periment, we obtain Ez ~ 176 K at a moderately strong 
field B = 30 T. Therefore, a typical value of the excitonic 
gap is A ~ 657 K. 

It was assumed in I that the gap A in Eq. (p~2|) is driven 
by the singlet excitonic order parameter 

(*<7 * } = ^KAa^KAa + ^ K , A ^K>A* 

a=±l 

-^KB^KBa-^ KIBa ^K'BcA), (13) 

where we explicitly wrote unit matrix oo acting on spin 
indices to underline the singlet character of this or- 
der parameter. The Zeeman interaction, however, may 
favo r 25 ' 31 the triplet order parameter, (^>a^), where the 
vector a = (ci, 02, 03) made from Pauli matrices acting 
on spin indices. Therefore, in addition to the singlet or- 
der parameter (fT3|) . one may also include the triplet order 
parameter— 

(^crs*) = <j{[^ KAa ^KAa + ^K, Aa ^K'A,j 

<T=±l 

-IpKB^KBa ~ IpK'Ba^K'Ba])- (14) 

In the language of symmetry, the Zeeman term explicitly 
breaks the Z7(4) symmetry down to the U(2) c x U{2)d (see 
Appendix C of J). [Note that the condensates (^fai^f) 
and (\&<72 ^) would spontaneously break the SO(2) sym- 
metry of the in-plane rotations and will not be considered 
here.] The dynamical generation of the gaps connected 
with order parameters (| 13[) and (fT4|) leads to the sponta- 
neous breakdown of the U(2) c x U(2)d symmetry down 
to the abelian U(l)i x U(l) 2 x U(l) 3 x E/(l) 4 one. There 
might exist also other gaps that break the U(2) c x U(2)d 
symmetry. Therefore, the problem is very complicated 
in general due to the potential possibility of many com- 
peting order parameters. 



III. PHENOMENOLOGICAL MODEL 

There are several ways to tackle the problem of com- 
peting order parameters in graphene. The best way 
would be to include all relevant interactions and study 
the instability of the system with respect to the formation 
of all possible condensates. An approach of this type was 
recently attempted in Ref. [25| in a weak-coupling regime. 
The applicability of such an approach is extremely lim- 
ited at strong coupling, however. 
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The second option is to use a phenomenological expres- 
sion for the Zeeman splitting as an input for the ther- 
modynamic potential and minimize the latter only with 
respect to the gaps of interest. In I we used the latter 
approach, taking as an input the usual Zeeman splitting 
([3|) and considering the thermodynamic potential with 
the singlet excitonic gap A in the weak-coupling regime 
with g ~ 0.07. 

The analysis is reliable in this regime, but the results 
do not describe all the features of the Hall conductiv- 
ity in graphene quantitatively. In fact, the analysis of 
the data presented in Sec. [Til led us to the conclusion 
that the weak-coupling regime is most likely improbable. 
Therefore, here we insist on the strong-coupling regime 
to describe the physics of graphene. Our argument is 
the following. The single t gap A in I was related to the 
Landau scale L(B) — \Jhv 2 F \eB\/c via the solution of 
the gap equation, A = bL{B), where the dimensionlcss 
parameter b in the LLL approximation is given by 
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V2 



dk e 



1 + k X o ' 



(15) 



with xo — 0.56y/2ng. Since the estimate for the value of 
A obtained in Sec. [TTJ is by an order of magnitude larger 
than in I, the corresponding parameter b = 0.4 is also 10 
times larger than in I. Then, using Eq. (fT5)l , we estimate 
that the value of the coupling constant is g ~ 1.56 which 
implies a strong-coupling regime. 

In this regime, unfortunately, there exist no reliable 
schemes for treating the pairing dynamics quantitatively. 
As was emphasized in Ref. HH, the problem is due to non- 
decoupling of the low-energy dynamics on the LLL from 
the dynamics on higher Landau levels. 

A general insight into the dynamics of magnetic catal- 
ysis at strong-coupling can be gained from the original 
work>2i in models with short range (contact) interactions 
in the leading order of 1/N expansion, where N is the 



number of fermion "flavors" (N equals 2 in graphene). 
The dynamical picture following from that analysis is the 
following. While an external magnetic held strongly en- 
hances the pairing dynamics, a nonzero density of carri- 
ers tends to suppress pairing. Moreover, for a wide range 
of parameters, the gap closes at the critical density that 
corresponds to filling the LLL. In terms of the chemical 
potential, this gives the critical value, \(i c \ — A in the 
absence of the Zeeman splitting^ and |/z c | = A + Ez 
when the Zeeman splitting is relevant. 

In addition to the exciton condensation in graphene, 
one should also account for the anomalously large Zee- 
man splitting (jlip . which is likely to have a similar dy- 
namical origin. As is discussed in I and in Sec. [TT] above, 
it is the only way to remove completely the four-fold de- 
generacy of the lowest Landau level necessary for the ex- 
planation of the experimentally observed v = 1 plateau 
in the Hall conductivity. 

Here, therefore, we implement all the details of the 
fermion pairing dynamics by making use of a simple phe- 
nomenological approach. In particular, we will use the 
simplest possible ansatz for the dynamically generated 
gaps A± for the spin up and down states, which cap- 
ture the essential features of the magnetic catalysis phe- 
nomenon discussed above. We write 



A_( J B, M ) = A(S,/i_), 



(16) 



where, as already defined in the text after Eq. (2}, 
fi± = fi ± Ez(B), with the Zeeman energy Ez(B) given 
in Eq. pip . Note that the triplet channel is taken into 
account by considering the two different gaps A± cor- 
responding to the order parameters 1/2(v1>(cto ± Vi)^), 
respectively. The magnetic catalysis dynamics is mainly 
realized at the two Fermi surfaces, fi± — 0, connected 
with the spin up and down quasiparticles, and leads to 
the gaps of approximately the following form: 



A(B,/x±) 



An 



0(B - B c 



B 

~B~r. 



arctan 



1 



arctan 



(17) 



Here B c is the critical field, 7 is the LLL width (or quasi- 
particle scattering rate) , and the dependence of the crit- 
ical chemical potential on the field B is 



fi c (B) = A 9(B - B c ) 



B 

K 



1. 



(18) 



The expression (TT7)) incorporates the following key fea- 
tures of the magnetic catalysis dynamics: 

(i) In view of expression ijTS]) . the gap is negligible for 
chemical potentials larger than the critical value 



fJ. c (B) which corresponds to the fil ling of the LLL, 
i.e., [i c ~ A(B, ti)]^ (see Ref . ■ 

(ii) In accordance with previous studies, 29 the magnetic 
catalysis occurs only when the field exceeds a crit- 
ical value B c . 

(iii) The effects due to non-vanishing scattering rate 7 
are incorporated through the smearing of the crit- 
ical region around fi c . Note that the choice of the 
arctan-function in the /u-dependence in Eq. (|17p is 
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suggested by the T = gap equation (All) in I. 

In total, we have three parameters, Ao, B Cl and 7 for 
the description of the gap generation, which can be used 
to fit the experimental data. In order to compare the 
theory and the experiment, we substitute the gaps A± 
given in (flB")) [with A(B,[i±) from Eq. (fT?) ] and depen- 
dences fJ,(V g ) and E Z (B) in Eqs. (JU) and (fTTj) (extracted 
from the experimental data) into the expression for the 



Hall conductivity: 

<J X y = ha xy (A + ,n + ) +a- av (A_,/x_)]. (19) 

The conductivity a xy in the limit B — > 00 can be 
expressed in terms of the digamma function ^(x) = 
d\nT(x)/dx (see Appendix B in J), namely 



o- xy (A±,n±) 



2e 2 sgn(eS) 



InW <5 



7tl . + i(ij.± + A ± ) 1 
2-kT 2 



It 



2ttT 



7tr+»(^±+A ± ) 1 

2ttT 2 



(A± 



r 



-A±)|, 
(201 



where 7t r is the transport scattering rate. One can verify 
that for T = 7tr = and A_ = A+ = A, Eq. (fT9|) 
reduces to Eq. @. 



IV. RESULTS AND THEIR INTERPRETATION 

The best fit is obtained while using the following values 
of the parameters: 

A = 680 K, (21a) 
B c = 7 T, (21b) 
7 = 7tr = 50 K. (21c) 

Note that these parameters capture essentially all the 
non-perturbative physics of a strongly interacting model 
at hand. Remarkably, it is sufficient to get a nearly per- 
fect fit for a wide range of the magnetic fields, see Fig. [4] 
In the upper panel of the Fig. [3J we present the experi- 
mental data for a xy , and in the lower panel their descrip- 
tion within our phenomcnological model. 

The theoretical curves plotted in Fig. [4] at T = 30 mK 
remain practically the same over a very wide range of 
temperatures up to T < Ao when 7 and 7 t r do not de- 
pend on T. At present there are no available experimen- 
tal data for the temperature dependence of v — 0, ±1 and 
v = ±4 plateaus. Furthermore, to compare our model 
with such experimental data at high temperatures one 
should also take into account the dependence of 7 and 
7tr on the temperature. 

The large value of Ao in Eq. (I21a|) . which is necessary 
to fit the data, and the simple estimate of the coupling 
constant g made below Eq. (|15p clearly imply a strong- 
coupling regime. The coupling constant is g = e 2 /cHvf = 
ac/eVF, where a ~ 1/137 is the fine-structure constant. 
With air on one side of the graphene plane and SiC>2 
on the other, the unscreened dielectric constant of the 
medium is estimated in Ref. Hi] to be e w 1.6eo. This 
corresponds to g w 1.37. We see that once again this 
supports the arguments that a strong-coupling regime 
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FIG. 4: (Color online) Hall conductivity from the data in 
Ref. |3 (upper panel) and in the theoretical model (lower 
panel) for magnetic fields B = 9 T (circle), 11.5 T (pentagon), 
17.5 T (hexagon), 25 T (square), 30 T (diamond), 37 T (up 
triangle), 42 T (down triangle), and 45 T (star). The pa- 
rameters in the model are 7 = 7 ti = 50 K, and temperature 
T = 30 mK. 



considered in this paper is more plausible than a weak- 
coupling one. 

It is instructive to consider how the quasiparticle en- 
ergy spectrum is affected by the applied gate voltage 
in graphene. Such a spectrum for a fixed value of the 
magnetic field is shown in Fig. [S] In essence, it is the 
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FIG. 5: (Color online) The theoretically reconstructed spec- 
trum of one-particle excitations in graphene at a fixed value 
of the magnetic field, B — 30 T, as a function of the gate 
voltage. The energy levels are shown by solid lines of differ- 
ent thickness which represents the degeneracy of the Landau 
levels: thin and thick lines denote single and double degen- 
eracy, respectively. The dotted line shows the dependence of 
the chemical potential on the gate voltage, see Eq. The 
points of the intersection of the chemical potential with the 
energy levels marks the position of the steps between neigh- 
boring plateaus in the corresponding dependence of the Hall 
conductivity on the voltage. 



of the v — for small /x (or V g ~ Vq) is the excitonic gap 
A, while the v — ±1 plateaus are related to the Zeeman 
term (see the discussion after Eq. (fT2")) above and Fig. 1 
in T). In the present scenario with a strong-coupling dy- 
namics, there are two gaps, A±, and for T = j tI = 
Eq. (fTT)|) reduces to 

e 2 

a xy = --sgn(eS) ^ 
x [sgn( M+ )0(|/x + | - A+) + sgn( M _)0(|/i_| - A_)] 

instead of a more simple Eq. ([2]). For the case Ez < A± 
shown in Fig. [5l this dependence of a xy on fx implies 
that while the size of the Oth plateau is <5/x = |A + + 
A_ — 2Ez\, the size of the ±1 plateaus is 8li\ = |A_ — 
A + + 2Ez\- When the effects of non- vanishing widths 7 
and 7t r are taken into account, the results become more 
complicated of course. However, the main qualitative 
features regarding the v = 0,±1 plateaus are captured 
already by Eq. ([22]) . It is clear, for example, that these 
plateaus result from a subtle interplay between the gaps 
A± and the Zeeman energy. We also note that the v = 
±1 plateaus disappear when A± = 0. It is also absent 
when Ez = and A + = A_ (see also the next section). 



strong-coupling pairing dynamics that is responsible for 
the drastic change of the energy spectrum at low volt- 
age. The four- fold degeneracy of the LL level is lifted at 
low voltage because of the gap formation and the large 
Zeeman splitting. When the voltage (and, therefore, the 
chemical potential and the density) is large, the pairing 
plays no important role. 

At large |V^|, the LLL, which has the energy Eq = 
and the four-fold degeneracy in absence of Zeeman split- 
ting, splits into two levels with the energies E = ±E Z 
and with the two- fold degeneracy (thick lines), as shown 
in Fig. [5] When the absolute value of /i is such that 
/x_ = n — Ez is smaller than the critical value fi c (B), 
the gap A_ opens. This gap causes the level with 
E = Ez to split into two nondegenerate levels with ener- 
gies E = Ez ± A_. Similar splitting of the energy level 
E = —Ez into two levels with energies E = —Ez ± A + 
is caused by a nonzero A + when the value = fx + Ez 
becomes smaller than ijl c (B). As should be clear from 
Fig. [5l it is only the presence of such nondegenerate lev- 
els in the energy spectrum of graphene that make the 
observation of the v = ±1 plateaus possible. In ap- 
proximately the whole region where both gaps A± are 
nonzero, the Hall conductivity of graphene develops the 
v = plateau. Moreover, this correspondence might be 
exact if the phase transition with respect to V g (or fi) is 
a strong first order phase transition. 

It should be noted that the mechanism behind the cre- 
ation of the v = 0, ±1 plateaus in the strong-coupling 
regime is different from that at weak coupling. In the lat- 
ter case, the dominant term responsible for the creation 



V. HALL CONDUCTIVITY WITHOUT 
ENHANCED ZEEMAN SPLITTING 

In order to appreciate the role of the enhanced Zeeman 
splitting, which is required to fit the data in Ref. 0, it will 
be instructive to consider also the case without Zeeman 
splitting, i.e., Ez — 0. It should be clear that this ap- 
proximation should also describe well a system with an 
unenhanced Zeeman energy Ez, as given in Eq. ([3]), when 
Ez < 7- By making use of the same model parameters 
as in the previous section, see Eq. (f2"Tj) , we plot the the- 
oretical curve for the Hall conductivity as a function of 
the gate voltage in Fig. [SI 

To calculate the dependence of the Hall conductiv- 
ity on the gate voltage, here we used the expressions 
in Eqs. (3.14) and (3.15) in Ref. [l4j. The correspond- 
ing result includes the contributions of all Landau levels, 
which are needed in order to reproduce not only the low- 
est plateaus with \v\ < 2, as in Fig.[U but also the higher 
ones. 

As we see from Fig. [6j the plateaus with v = ±1 and 
v = ±Ak (with integer k) do not appear when the Zee- 
man energy is negligible (or Ez ^7). Of course, this 
conclusion is hardly surprising in view of the correspon- 
dence between the energy spectrum in graphene and the 
form of the Hall conductivity illustrated in Fig. [5] When 
the Zeeman splitting of the Landau levels disappears, the 
plateaus with v = ±1, v — ±4, v — ±8, etc. collapse 
into a point. In the language of symmetry, this picture 
corresponds to a partial removing the degeneracy of the 
LLL, when the U(4) symmetry is broken down to the 
U(2) a x U(2) b (see Appendix C in J). 



9 



A =680 K y=50K 
~ B C =7T y, r =50K 


n^f — J n 


mm 

\ li 


1 ffmf 






1 



-75 -50 -25 25 50 75 
V g [V] 



FIG. 6: (Color online) Hall conductivity in the case of a negli- 
gible Zeeman energy. As in Fig. [4] the symbols label different 
values of the magnetic field: B — 9 T (circle), If .5 T (pen- 
tagon), f7.5 T (hexagon), 25 T (square), 30 T (diamond), 
37 T (up triangle), 42 T (down triangle), and 45 T (star). 



VI. CONCLUSION 

In this paper, based on the experimental data of Ref.H, 
we developed a model for the exciton condensation dy- 
namics in graphene in a strong magnetic field. The non- 
perturbative dynamics corresponds to a strong coupling 
regime with a realistic value of the coupling constant 
g — e 2 /ehvp of the Coulomb interaction. 

On theoretical side, the model incorporates the main 
features of the phenomenon of the magnetic catalysis-'' 
(the generation of an excitonic gap in a strong mag- 
netic field) and its realization in graphene considered 
in Refs. |28||29| . On phenomenological side, the model is 
based on features extracted from the experimental data 
in Ref. 0. Among them, the relation between the gate 
voltage applied to the graphene device and the chemical 
potential fj, is particularly important. Another important 
point established is a strong enhancement of the Zeeman 
splitting (spin gap) in graphene in strong magnetic fields 
B > 9 T. The existing experimental data, however, do 
not allow to extract the dependence of the spin gap on 
T and /i and here we considered its dependence on the 
field B only. 



The enhanced Zeeman splitting alone does not allow 
to explain the occurrence of the additional plateaus with 
the filling factors v = ±1 in magnetic fields B > 20 T ob- 
served in Ref. 0- In addition to such a splitting, it is nec- 
essary to remove the sublattice degeneracy in graphene 
to explain the origin of these plateaus. The available 
experimental data£ already contain a lot of constraints 
on possible microscopical mechanism which removes this 
degeneracy. In addition to this, the observation of the 
v = plateau at rather high temperatures^ indicates 
that a strong coupling regime is more believable than 
the weak coupling regime. 

Still further experiments are necessary to establish the 
temperature evolution of the v — 0, ±1, ±4 plateaus and 
their sensitivity to the quality of the samples. This will 
hopefully allow to decide which of the mechanisms men- 
tioned in the Introduction is realized in graphene in a 
strong magnetic field. For example, there is a suggestion 
how to detect a gap either in microwave^ or optical 3 -^ 
response. Recent measurements done in the far infrared 
region 38,39 show that the second method may work if the 
measurements will be done not in epitaxial graphite as 
in Ref. [38| or in highly oriented pyrolytic graphite as in 
Ref. [39] but in graphene in a strong magnetic field. 
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